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ABSTRACT 

We present a method to characterize the spectral transfers of magnetic energy between scales in 
simulations of stellar convective dynamos. The full triadic transfer functions are computed thanks to 
analytical coupling relations of spherical harmonics based on the Clebsch-Gordan coefficients. The 
method is applied to mean field ad dynamo models as benchmark tests. From the physical standpoint, 
the decomposition of the dynamo field into primary and secondary dynamo families proves very 
instructive in the aQ case. The same method is then applied to a fully turbulent dynamo in a solar 
convection zone, modeled with the 3D MHD ASH code. The initial growth of the magnetic energy 
spectrum is shown to be non-local. It mainly reproduces the kinetic energy spectrum of convection 
at intermediate scales. During the saturation phase, two kinds of direct magnetic energy cascades are 
observed in regions encompassing the smallest scales involved in the simulation. The first cascade is 
obtained through the shearing of magnetic field by the large scale differential rotation that effectively 
cascades magnetic energy. The second is a generalized cascade that involves a range of local magnetic 
and velocity scales. Non-local transfers appear to be significant, such that the net transfers cannot 
be reduced to the dynamics of a small set of modes. The saturation of the large scale axisymmetric 
dipole and quadrupole are detailed. In particular, the dipole is saturated by a non-local interaction 
involving the most energetic scale of the magnetic energy spectrum, which points out the importance 
of the magnetic Prandtl number for large-scale dynamos. 

Subject headings: Dynamo — Magnetohydrodynamics (MHD) — Stars: magnetic field — Turbulence 



1. INTRODUCTION 

Magnetic fields are observed in astrophysical bodies 
in a broad range of scales, from the object scale to th e 
smallest dissipative scales (Donati & Landstreet 2009). 
The origin of such fields is, in most cases, due to a 
hydromagnetic dynamo process. Recent developments 
in dynamo theory led to a disti nction between large- 
scale and small-scale dynamos ( |Cattaneo fc Hughes | 
2001|). The large-scale dynamos produce magnetic 



fields at larger scale than the largest velocity scale (or, 
the largest driving scale), while small-scale dynamos 
generate magnetic fiel ds at all scales sm aller or equal 
to the driving scales ( |Tobias et al. 2011). Large scale 
dynamos also sometimes refer to dynamos that develop 
a large-scale magnetic fiel d in super-equipart ition with 
large scale kinetic energy (Olson et al. 1999). In that 
case, small-scale dynamos refer to those that develop 
a spectrum peaked at small scales. Both dynamos are 
generally acting together, like in the Sun, where we 
observe both large-scale, intense, global magnetic field s 



dSchrijver fc DeRosa|[2003l peRosji et al.||20lT^ 
and small-scale mag netic fields ([Hagenaar et al 



Centeno et al. 2007). In the case where multiple 
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scales 

coexist also in the velocity field, special care i s needed 
to isolate 'small ' and 'large' scales tendencies (Tobias & 
|Cattaneo||2008ap . 
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Recent developments in numerical simulations in 3D 
spherical geometry allow us to model fully non-linear 
dynamos in stars, involving a broad range of scales 
(|Brun et all |2004l |Browning| [20081 |Br own et al.||2010" 
R" 



acine et al.||20lT| [Kapyla et al. 2012). The now scales 
ol a stellar convection zone extend from the large-scale 
differential rotation down to the smallest convective 
scales. In order to properly characterize such dynamos, 
one may use specific methods to tackle the multi-scale 
aspect of the problem. The principal tools that have 
been used in the literature are spectral decomposition 
([Frick fc Sok oloff 1998] |Dar et al.|[2001[ ), and wavelet 
analysis (Farge 1992). In this paper, we choose to use 
spherical harmonics decomposition (which is adapted 
to the spherical geometry of stars, |Bullard fc Gel lman 
(1954)) to develop a spectral analysis of energy transfers 
in the frame of dynamo theory. 



Mainly used to st udy turbulenc e (|Frisch 1995 



Debli quy et~aT1 [2005] |Lesieur| [20081 |Alexakis etlif 
2005), spectral analysis is also a useful tool to char- 



acterize magnet o hydro dynamic (MHD) processes like 
dynamos (|Biskamp | [ 19931 [Blackman & Brandenburg 
2002[ |lVlminni et al.lTM75[^ivermore et al.||2010|) or 
the magne t o-rotat ional instability (MKl, see |Lesur fc 
Longaretti ( 2011[ )). Understanding spectral energy 
transfers between scales in such processes may reinforce 
our ability to characterize non-linear MHD phenomena. 
The shell-to-shell or mode-to-mode methods have been 



2 



Strugarek, et al. 



recently and extensively used in the context of MHD 
turbulence. Indeed, the classical Kolmogorov approach 
to turbulence must be adapted to the MHD case, since 
the magnetic field in duces an anisotropy that has to be 
taken into account fllroshnikovl 119641 iKraichnanl 11965 



Biskamp|)1993| |Goldreich jF^mEar]|1995p . BepencT- 
ing on tne dimensionality of the problem, spectrum 
slopes are often understood to result from local (direct 
or indirect) transfers of energy, r eferred as cascades 
(|Biskamp||l993| |Maron et~aT1 |2QQ4j ). However, it was 
found that non-local interactions m MHD turbulence 
may also contribute importantly to the built-up of th e 
spectrum ( [Schilling fc Zhou||2QQ2| [Aluie fc Eyink||2010| ). 
The directions and localizations oi energy transfers are 
then less obvious to identify, and studies dedicated to 
transfer processes in spectral space are essen tial to prop- 
erly understand spectrum slopes in MHD dPolitano fc 
Pouquet|199"8l|Boldyrev et al.|2009||Pouquet et al.|2Ulip . 



In the past, spectral analyses have mainly been 
used with Fourier spectral decomposition, generally 
in cartesian coordinates and periodic parallelepipedic 
boxes. The Fourier decomposition is indeed a natural 
way to understand spectra, since the Fourier wave 
numbers represent the inve r se of a spatial scale. More 
recently, Hughes & Proctor (2012) used Fourier spectral 
analysis to study the influence of large-scale sheared 
flows on local convective dynamos. They show that 
the dynamo process depends on a broad range of scales 
in this case. We also study dynamo process in the 
present paper, though in spherical geometry with a 
self-consistently generated large scale sheared flow (the 
differential rotation). In the case of stars or planets, 
the spherical geometry of the object makes the spherical 
harmonic s basis much more ada pted to the spe ctral 



analysis (pullard fc Gellman 



1954[ ). For example, livers 



& Phillips (2008) wrote the decomposition of the MHD 
equations onto spherical harmonics in the framework of 
geodynamics. They were able to analytically express the 
non linear terms by calculating the coupling between 
the spherical harm onics with the Cle bsch-Gordan 
coefficients (see also [Mathis & Zahn 20051 in a stellar 

. \ rrvi l I 1— I ■ ' 



context). The spherical harmon i cs aec omposition was 



also used by 



(2010) to develop a 



Liver more et al. 
spectral analysis which is similar to the one we present 
in this paper. They used it to identify the spectral 
interactions leading to a different saturation level of 
large-scale magnetic field in kinematic and non-linear 
forced dynamos. In the latter case, they observe a 
significant reorganization of the magnetic field such 
that a strong large-scale magnetic field can emerge. In 
addition to a particular geometry, the choice of a certain 
basis for the spectral analysis may be motivated by the 
presence of anisotropy (e.g., between the vertical and 
horizontal directions), which is poorly de scribed by the 



classical Fourier decomposition (e.g., see Rincon 2006 
in the case of turbulent convection). 

The spectral interactions in MHD involve triads 
coupling, meaning that two modes interact to impact 
a third one through a triangulation rule. Depending 
on the ideal MHD invariant considered, these kinds of 
interaction involve couplings between the velocity and 
the magnetic field, impacting the magnetic or the veloc- 



ity field. Shell-to-shell methods generally only consider 
dual interactions, raising an ambiguity on the medium 
(third com ponent) of the triadic interaction (Verma 
et al.|[2005| . In order to cope with thi s ambiguity, other 
studies (e.g., [Schilling fc Zhou||2002[ ) made use of the 
eddy-damped quasi- neutr al Markovian (EDQNM) two 
point sta tistical closure (Frisch et al. |1975 Pouquet 
et al.||1976[ ) to get an analytical expression of the triadic 
interactions. We point out that such methods are rele- 
vant for, e.g., developing subgrid-scale models f or large 
eddy simulations (LES). |Tobias fc Cattaneo ( |2008b[ ) 
demonstrated that such truncated methods often badly 
describe coherent turbulent structures in flows, which 
are thought to be responsible for the generation of 
large-scale fields in dynamos. This limitation is relevant 
for reduced spectral models, which aim to reproduce the 
full turbulent behavior with a reduced number of modes. 
However, here we directly calculate the full triadic 
shell-to-shell interactions of all the scales included in 
our simulations, (i.e., we do not use any specific closure 
in spectral space to compute the full triadic interactions). 

We here applied our spherical harmonic based method 
to a numerical simulation of a stellar convection by con- 
sidering the three spectral components of the triadic in- 
teractions. The originality of the method we develop 
in the present work resides in the facts that (i) we de- 
compose explicitly the spectral interactions for both the 
magnetic and velocity fields, (ii) calculate explicitly all 
the coupling coefficients between those fields and (Hi) we 
use it to study dynamo action in a solar-like turbulent 
convection zone that possesses self-consistent large-scale 
flows (differential rotation, meridional circulation, ...) as 
well as a broad range of turbulent scales. 

In section [2] we present the set of MHD equations we 
will use, derive from them the spectral evolution equa- 
tion for the magnetic energy in the spherical harmonics 
formalism and analytically validate our method. A toy 
model of an axisymmetric a — ft dynamo is analyzed with 
our spectral method in Sect. [3] In section [4j we apply 
our method to study non-linear dynamo action in a nu- 
merical simulation of a solar convective zone. Finally, 
conclusions and perspectives are given in section [5] 



2. MAGNETIC ENERGY EVOLUTION EQUATION 

2.1. Main equations in physical space 

We use the well-tested Anelastic Spherical Harmonics 
(ASH) code which models turbulent stel lar convection 
zones ( |Clune et al.|1999 Jones et al.||2011 ). It solves the 
following three di mensional MHD set of equations (see 
|Brun et aL ( 2004[ )) in the anelastic approximation, in a 
reference frame rotating at the angular velocity fto = 



Spherical spectral transfers: magnetic energy 



3 



ftoe z (where e z is the cartesian vertical axis): 

V • (pU) = 0, 
V-B = 0, 

p [dtU +(U-V)U + 2O xU] = - VP + pg 
1 



(1) 
(2) 



47T 



(V x B) x B 



pf [d t s - 



U- V(5 + 5)] : 
+K pT VS - 

+2pp[eijeij 



V • D — [VP — pg\ , 

(3) 

V • [ Kr pc p V (T + T) 



npTVS] 

1 



(v-un. 



d t B = V x (U x B) - V x (t?V x B) . 



(4) 
(5) 



where the spherically symmetric background thermody- 
namical state is denoted by bars (fluctuations with re- 
spect to the background state are denoted without bars), 
v is the local velocity, Kjf. IS the radiative diffusivity, and 
ft, v and r] are respectively the effective thermal diffu- 
sivity, the eddy viscosity and the magnetic diffusivity. 
The thermal diffusion coefficient applies at the top of 
the modeled convective zone (where convective motions 
vanish), to ensure the heat transport through the upper 
surface. J = (c/47r)V x B is the current density, and the 
viscous stress tensor D is defined by 



Dij = 



-2pv 



1 



(V • U) 6 i: 



(6) 



where e^- is the strain rate tensor, and Sij is the Kro- 
necker symbol. The system is closed by using the lin- 
earized ideal gas law: 



P 
P 



T 
T 



P 

IP 



S_ 



(7) 



with c p the specific heat at constant pressure and 7 the 
adiabatic exponent. The vectorial fields are decomposed 
in poloidal and toroidal components: 

B(r, ^) = VxVx [C(r, 0, ip)e r ] 

+ Vx [A(r,6,<p)e r ], (8) 
p(r)U(r,0,y>) = V x V x [W(r,0,ip)e r ] 

+ V x [Z(r,e,(p)er] , (9) 

where (e r ,e<9,e^) are the unit vectors in spherical co- 
ordinates. All the quantities are time-dependent. This 
decomposition ensures numerically that both the mag- 
netic field and the mass flux remain divergenceless up to 
the machine precision. 

A potential match of the magnetic field (V x B = 0) is 
applied both at the bottom and top radial boundaries. 
For the convective dynamo case (Sect. [4|, the boundary 
conditions for the velocity are impenetrable and stress- 
free. A latitu dinal entropy gradie nt is imposed at the 
bottom (as in |Miesch et al. (2006)), and we fix a con- 
stant entropy gradient at the top of the domain. 

2.2. Magnetic energy transfer functions 

2.2.1. The formalism 

In this section, we present a method to obtain a spec- 
tral (in the sense of the spherical harmonics) evolution 



equation for the magnetic energy, starting from the in- 
duction equation ([5f . In order to deal w ith vectorial 
fields and spherical harmonics (see Eq. (A2) in Ap- 
pendix), it is pr actical to define the vectorial spheric al 



pendixj, it is practical to denn e tne vectorial spneric a 
harmonics basis (Rieutord 1987; Mathis fc Zahn||2005 ): 



■(0,<P) 



V ± Y t 
V ± 



eg 



x R 



-J—f) V 

sin9 U( P 1 l 



r) Y 71 



(10) 

It is an orthogonal basis for the scalar product JL • dfi, 
where S is a spherical surface and d^ = sinOdvdip the 
associated infinitesimal solid angle. The mode numbers 
m and I are the azimuthal wave number and the spherical 
harmonic degree (which characterize to their latitudinal 
variations). The g enera l properties of this basis maybe 
found in appendix |A.1| The two main vectorial fields 
that appear in the induction equation ([5| are the mag- 
netic and the velocity fields. We want to project those 
fields on the vectorial basis (10), using the decomposi- 
tions ([8|-([9|. Fortunately, the curl of a vector is a linear 
operation that can be expressed very easily i n th e vec- 
torial spherical harmonics basis (see equation (A9)). We 
obtain: 



1=1 m=-l k 



B(r,e,<p) 



d r C l m (r)ST+ ^^tA, (11) 



p(r)U(r,e,f) = £ E |^^<(r)RF 

1 = 1 m = -l 



9 r Wi(r)S™ + ^^Tf 1 1 , (12) 



where we have projected the toroidal and poloidal com- 
ponents of the fields (A, C, Z and W) on the scalar 



spherical harmonics basis. We see that in equation (11) 



the poloidal C and toroidal A components of B are re- 
spectively projected on (R, S) and T. Consequently, in 
the remainder of this paper the projection of any vecto- 
rial field on (R, S) will be referred as poloidal, and the 
projection on T as toroidal. 

2.2.2. Shell to shell analysis 

To study the transfers of energy between scales on a 
spherical surface, we distinguish the different scales of the 
axisymmetric (m = 0) and non-axisymmetric physical 
fields by defining shells L° and L* as follows: 



■C l m T?}. 



(13) 
(14) 



-l<m<l 



This distinction is natural when studying the generation 
of large scale axisymmetric field. Another choice of shells 
based on dyna mo f amilies will also be used in this paper 
(see Appendix A. 7 and the end of Sect. 2.2.3). Note that 
the defined shells are orthogonal, i.e. that any scalar 
product of strictly different shells is zero. In order to 
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simplify the notations, the shell L may represent either 
axisymmetric or non-axisymmetric shells. Exponents 
and * denote axisymmetric and non-axisymmetric com- 
ponents in the remainder of this paper. 
The shells of magnetic energy in spectral space are then 
defined by 



This operation is not easily calculate d in the (R, S,T 



lis 



Br- Br dfi. 



(15) 



This spectrum may be defined at any radial location of 
the spherical domain (e.g., in a stellar interior). Our 
spectral analysis intends to characterize horizontal scales 
of velocity and magnetic fields, and can be easily applied 
at any and as many as necessary depths, depending on 
the radial locations one wants to focus on. In addition, 
the different terms of the evolution equation of the mag- 
netic energy (see next section) explicitly depend upon 
the radial gradients of the different quantities. The hor- 
izontal couplings produced by the vertical interactions 
are thus taken into account by our description. Finally, 
it is worth noting that the spectra may be different at 
various depths in a convective dynamo model. In this 
work, we will only study pure convection zone dynamos, 
we will consequently focus on the spectral interaction in 
the middle of the convection zone. 

2.2.3. Spectral magnetic energy equation 

In order to obtain the spectral magnetic energy evo- 
lution equation, we multiply equation (5|) by and 
integrate it over the spherical surface so 



lat 



dtET* 



Vj 



{V L (L U L 2 )+F L (L U L 2 )} , (16) 



where V regroups the diffusion terms, V represents the 
volumetric production of magnetic energy, and T is the 
divergence of the flux of magnetic energy through a 
spherical surface. Note that the production V term has 
to be understood as a general production term, that can 
either be positive (real production) or negative (destruc- 
tion). The sum over L\,L 2 involves the triangular se- 
lection rule — l 2 \ < I < h + fa, ttii + m 2 = m} that 
comes natural ly fro m the spherical harmonics coupling 
(see Appendix A.2). The expressions of the three contri- 
butions X>, V, and J 7 are given by: 



(17) 
dft, (18) 



%(r)= / {r?B L -V 2 
Js 

+d r r]e r ■ (B L x V x B L )} dfi, 

V L (r,L u L 2 ) = [ (U Ll xB L2 )-VxB L 
Js 

T L (rMM)= [ V-[(U Ll xB L2 )xB L ] dfi. (19) 
Js 

We have split the diffusive term by considering a diffu- 
sivity r] that only depends on r. The production and flux 
of magnetic energy are discretized over scales so that we 
compute which scale of the velocity field (L\) is interact- 
ing with which scale of the magnetic field (L 2 ) towards 
a studied scale L. 

Although the expressions are formally writ- 

ten, they include the evaluation or vectorial products 
decomposed on the vectorial spherical harmonics basis. 



basis thus we use an alternative basis ( Varshalovich et al. 
119881) to compute it. For sake of simplicity, these details 



are given in appendix |A. 5 

Alternatively to considering axisymmetric and non- 
axisymmetric spectra, it is instructive to decompose the 
flow and field into the so-called primary (dipolar, an- 
tisymme tric) and secondary (quadrupolar, symmetric ) 
families (Mc Fadden et al.||1991| [Roberts fc Stix||1972[ ). 
These families were proven very insightful t o character- 
ize geophysical and astrophysical dyn amos ( |Gubbins fc| 
|Zhang|1993l [DeRosa et al.|201 1[ [2012] ) . Further, the pri- 
mary/secondary distinction greatly simplifies the trans- 
fer maps of V and T . Indeed, the coupling between 
fields of the same family always gives secondary fields, 
while the coupling between fields of differ ent fa milies al- 
ways gives a primary field (see Appendix A. 7). Both of 



these approaches (i.e., axisymmetric/non- axisymmetric 
and primary/secondary distinctions) will be used here- 
after. 

2.3. Validation and illustration of the method 

In this section, we illustrate the coupling calculations 
for two simple fields. The reader only interested in physi- 
cal discussions may skip this part and go directly to Sect. 

® 

Since the ASH code is a spectral code, it solves the MHD 
equations for the spherical harmonics coefficients of the 
fields. Although it does not compute explicitly the de- 
composition on the vectorial spherical harmonics basis 
(10), it is straightforward to make use of this basis in 
the code by using the transformation relations ([TT])-( 12). 
We have added in the code the ability to compute the 
different terms (17)-(19) of the spectral magnetic energy 
equation (16). 



In orde r to ill ustrate and validate both the coupling co- 
efficients ( A17) for the vectorial product and the general 



method, we numerically computed a simple analytical 
test case. We initialize the magnetic and velocity fields 
in the following way: 

B = a (b? + isi - iltr 1 ) + » (s? + 5S! - isr 1 ) 



(i? + iTi-iTr'). 



pU = d (R4 - R-2 



(s 2 



S2" 1 ) 



7(12 



where a, 6, c, d, e and / are functions of r only. 
This initialization allows us to test at the same 
time the axisymmetric/non-axisymmetric and non- 
axisymmetric/non- axisymmetric coupling schemes be- 
tween the velocity and the magnetic fields. The low 
order harmonics (I G {0,1,2}) that are involved make 
the analytic calculation easy. We display on figure [l] the 
possible couplings (via vectorial product) between the U 
and B fields we initialized. This is in fact a schematic 
representation of the triangulation rule that appear in 
the summation of equation (A17). The analytical cal- 
culation of the values of the three large green circles is 
given in appendix|B] The resulting vectorial products cal- 
culated by the code using the Wigner coefficients show 
very good agreement with the coefficients calculated an- 
alytically (table [l] in appendix [B| . 
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Fig. 1. — Schematic coupling between spherical harmonics of U 
and B for the simple test case. The B modes are represented 
by blue squares, and the U modes by red triangles. The black 
arrows represent the coupling between the modes, the green circles 
on them represent resulting modes obtained from the coupling via 
the triangulation rule of the vectorial product U x B. The three 
highlighted green circles in the center correspond to the modes 
calculated in appendix [B] 



We stress here that this test has been done for low m 
and I values. The numerical accuracy of the algorithms 
calculating the Clebsch-Gordan coefficients (and thus 3j, 
6j and 9j Wigner coefficients) is known to decrease with 
increasing I and m. The calculation routines we use are 
accurate up to values of I of the order of 500. To do so, 
we used a multiple precision packag^Jto simulate large- 
precision numbers that are needed to compute the ratios 
of factorials and binomial coefficients that are involved 
in the Wigner coefficients calculations. However, the cal- 
culation time of the transfer functions V and T increases 
dramatically with I and m. For practical reason, when 
computing fully nonlinear dynamos (see Sect. [4] below), 
we have chosen to limit the computation of the coupling 
coefficients to / max = 70, even if the effective resolution of 
such simulations reaches / max = 340. From time to time, 
we do calculate the transfer terms for high Vs to have 
an indication of how energy is transfered at the smallest 
scales (see Sect. [4|. Nevertheless, the magnetic-energy- 
carrying scales in the spectrum are dominated by I < 70 
in this case. We thus capture the essential part of the 
dynamics. 

3. AXISYMMETRIC aft DYNAMO 

In this section, we use the spectral method we devel- 
oped in Sect. [2] on two aca demic cases. F irst, we explain 
how the classical ft effect (M offatt 1978) is represented 
by our formalism (Sect, |3.1| )~ Then, we calculat e th e 
spectral transfers for a mean field aft model (Sect. |3.2[). 



3.1. Omega effect 

The complexity of the two spherical harmonics bases 
may be confusing when it comes to interpret simple and 
classical dynamo processes. We thus give hereafter a 
step-by-step explanation of the ^-effect in the two vec- 
torial spherical harmonics bases formalism. 



We start with a purel y di polar poloidal magnetic field 
that reads (using Eq. (A9)) 



B p (r,0) = & r (r)R? + &*(r)S?. 



(20) 



Then, we want to calculate the effect of a differential 
rotation that reads 



ft = A + B cos 2 . 



(21) 



Such differential rotation is usually seen as a "Z = 2" 
field. Though, it projects on a I — 3 component when 
considering the azimut hal component o f the velocity 



- 

reads 



sin Oft (see Roberts & Stix (1972)), which 



U^(r,0) = ^(r) (A sin 
~U^r)d e Y 3 °e^-- 



B sin cos 9)e v (22) 



(23) 



In general, Eq. ( |22| ) should project both on and 
T3. For the sake of simplicity, we select here a pro- 
file of differential rotation that is purely described by a 
(I = 3 , m — 0) harmonic, which corresponds to B = 5A 
(Eq. pi))). We simply apply the cur l ope rator (A9) and 



make use of the coupling relations (A17) to obtain the 



production of B in the induction equation 



V x (Uy, x B p ) 



f(b r ,beMT 2 . 



(24) 



We recovered that the action of differential rotation on a 
purely axisymmetric poloidal field creates a toroidal field 
Bp oc sin#cos#. With our notations, this kind of field 
will be labeled as a '/ = 2' field. 

An additional feature of the differential rotation can 
also be learnt from this little analysis. We immediately 
remark that for axisymmetric fields, the first Wigner co- 
efficient involved in the coupling between t wo she lls L\ 
and L 2 is zero if 3 + l\ + I2 is odd (equation ( A18| )), i.e., 
if h and I2 are of the opposite parity. The shearing effect 
of differential rotation will then always couple axisym- 
metric scales of the magnetic field that are of opposite 
parity, which will be observed in the transfer maps in 
more complex cases (e.g., Figs[3[b) and [l3[ ). 

This simple example s trikin gly highlights how the vec- 
torial product formula ( |A17[ ) couples together two sim- 
ple fields. This d escri ption of the ft effect will guide our 
analysis in Sects. |3.2| and|4| 

3.2. Case of a cyclic mean field dynamo 



We use the ASH code flClune et alT]|1999[ |Brun et al. 
2004|) to simulate an axisymmetric mea n held dynamo 
( |Charbonneau||2010| |Jouve et aLl|2008| ). To do so, we 
solve only the induction equation considering uniquely 
a pre scribed differentia l rotation profile (see also section 
3.8 of |Jouve et al. 2008, for a similar use of a 3D spherical 
code to model a — U dynamos). 

Our radial domain is defined between r& = 0.6 Rq and 
r t = 0.966 R G . We use a resolution of N r x Nq x = 
64 x 128 x 256. We choose a solar differential rotation pro- 
file ftDR(r,0) = —doZjjR/^smO) through the toroidal 
component Zjjr of the momentum, which is, in the frame 
rotating at ft = 2.6 10~ 6 s _1 : 



1 http : / / crd- legacy . lbl . gov/~dhbailey/mpdist/ 



Z D R(r,9) = Z t (r) [ A cosd + - cos 3 



C 



(25) 



(a) 
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Br at r = 0.92 R 
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Fig. 2. — (a) Solar like differential rotation profile computed from |Schou et alT] ( |1998| >. The tachocline is located near r = 0.7 Rq, and 
the base of the tachocline rotates at the solar rotation rate f2o = 2.6 10 — b s~ 1 . (b) Butterfly diagram in the axisymmetric a — Q dynamo 
(only B r as a function of time and latitude is shown, in the upper convection zone). The vertical black dotted line represent the time at 
which we display the spectral interactions in Fig. [3] 



From Schou et al.| (|1998|) , we take A 
nHz and U 



329 



257 nHz, B = 321 
differential rotation then 



nHz. The 

naturally projects on Uj, U3 and U5. The radial profile 
Z t (r) is chosen such as to simulate a stable region at the 
base of the domain and is defined by 



Zt{r) = p- 



1 + tanh 



4.8710 



10 



210 9 



(26) 



We initialize our magnetic field with a seed I = 3 
poloidal (antisymmetric and axisymmetric) field. Finally, 
we add an a effect to the induction equation such that 



d t B = V x (U x B 



V x (77V x B) . (27) 



Since we do not take into account in this simple case 
the feedback of the Lorentz force on the flow via the 
Navier- Stokes equations, since we only solve the induc- 
tion equation, we need to quench the a effect. Hence, a 
is defined by 



a(r, 0) = ao e 



(- 



cost 



\B\/B q ) 



(28) 



This is the simplest a th at is needed to trigger an oscil- 



lating solar-like dynamo ( Charbonneau [2010 ) ; it is anti- 
symmetric with respect to the equator. The radial profile 
of a is localized near the base of the convection zone and 
the quenching value is given by B q = 10 3 G. We have 
deliberately chosen an a-effect that operates only on the 
poloidal component of the induction equat ion, therefor e 
computing an aft mean field dynamo flMoffattl [1978]) . 
This aQ dynamo exhibits the characteristic butterfly di- 
agram showed in Fig. [5|b) (at r = 0.92 Rq). Although 
this a profile is ad-hoc and one among t he many profiles 
that were tested in the literature (e.g., |Roberts & Stix 



19721 ICha rbonne au fc MacGr egor 1997; Bon a nno et a 
200^pnang et al.|2Q03|JJouve et al.[2008p , we chose this 
form because it easily triggers an oscillatory dynamo and 
its effect in spectral space can be easily calculated. It is 
consequently a good choice to illustrate our new spectral 
method. With the parameters we chose, the cycle period 
is of the order of 400 days (see Fig. |2jb)). 

The extra a effect adds a new term in the spectral 
energy equation (16) that can lead to complex formula 



in spectral space. We rewrite the energy equation 



/ V x (aB (p e (p ) L B L 
Js 

J2 {Pl(L 1 ,L 2 )+F l (L 1 ,L 2 )} • (29) 



The interested reader may read Appendix | A.4| for a com- 
plete spectral description of this a effect. 

Wherever |B| is not too large, the quenching part of 
the a effect is negligible. In that case, the a effect (which 
restores poloidal field from toroidal field) simply cou- 
ples a L shell of toroidal field to its neighboring shells 
of poloidal field, namely L — 1 and L + 1. When |B|/£? g 
becomes large, the a effect is quenched and the poloidal 
magnetic field stops being restored. When it is suffi- 
ciently low, a stops being quenched and the poloidal field 
grows again. This sets up a simple feedback mechanism 
and a cycle establishes. 

The magnetic energy spectrum is dominated by a 
L = 2 component that sets the phase of the total cy- 
cle. The various shells energy oscillate with roughly the 
same period, but are generally out of phase. This phase 
shift is a natural ingredi ent that allow the reversal of 
the o verall field polarity fKnobloch et al. 1998 Tobias 
^002] ). 

We stress here that the magnetic field created in this 
experiment is of the primary family ((R^+i? ^+1) and 
see appendix A.7). Our initial magnetic field is a 



poloidal primary field (R3, S3). As a result, the toroidal 
field created thr ough the ^-effect is also a primary field 



(T2, see section 3.1). Then, our a-effect, that creates 



poloidal field from toroidal field, transforms the primary 
toroidal field into a primary poloidal field (R^S^), 
which is of the same type than our initial magnetic field. 
Hence, no secondary field can be created in the simula- 
tion (which is confirmed by our results), and the a-effect 
can only act on the primary toroidal field to create a pri- 
mary poloidal field. This is a direct consequence of the 
well-known separability property of the induction equa- 
tion between the dipolar and quadrupolar families, when 
symmetric flows and ant isymmetric a effect are chosen 
dGubbins fc Zhang|[l9g3l ) . 
We display in Fig. |5Ja) the evolution of the different 



Spherical spectral transfers: magnetic energy 



7 




1000 1500 2000 2500 3000 2 4 6 8 10 

time [days] U axi 

Fig. 3. — (a) Evolution of the different terms of equation ( |29| for the shell L = 2. Production V is the dashed blue line, flux T the 
dashed red line, diffusion T>\ + T>2 the dotted black line and the a effect is the dash-dot green line. The black plain line is the total of all 
the contributions, (b) Flux T contribution to the L = 2 shell (see Eq. \29\ ) at t = 1500 days (vertical dotted line in panel (a)). The 2D 
color maps are the B° — U° transfer functions (dark red is the maximum value, black the minimum). The oblique black lines represent the 
boundaries of the triangular selection rule. 



terms of the magnetic energy equation ( 29 ) for L = 2 at 
r = 0.92 R Q during the same time perioa than the but- 
terfly diagram in Fig. f2^b). The primary toroidal field 
energy clearly evolves due to the production V (dashed 
blue line) and flux T (dashed red line) terms that account 
for the effect of differential rotation on the magnetic field 
(as expected, since there are no other production nor ad- 
vection terms). The two terms cancel each other out with 
a small time-lag, their sum combines with the ohmic dif- 
fusion (dotted line) to produce oscillations (solid line) of 
the total L = 2 energy. Note that the a effect plays no 
role et r = 0.92 Rq since it is concentrated at the base 
of the "convection zone" (equation ([28]) ) . 

Our new method allows us to characterize how scales 
interact to produce this behavior. We display in Fig. 
|3Fb) the transfer map for the flux F term of equation 
Tz9| for the L = 2 shell at its maximum. The differ- 
ential rotation is composed of the U?, U3 and U5 shells 
(eq. (25)). The transfer maps during minima (not shown 
here) are qualitatively opposite, which means that all the 
couplings between the shells reverse sign during the cy- 
cle. This reversal of all shells is a simple, direct conse- 
quence of the reversal of the whole magnetic field. At 
this position, the poloidal magnetic energy (not shown 
here) evolves because of the ohmic diffusion of the a- 
driven poloidal field at deeper radii. Here, the poloidal 
magnetic field couples with the differential rotation to 
transfer energy to the toroidal L = 2 magnetic shell. The 
U3 — appears to be the dominant interaction that sets 
the L = 2 cycle. Interestingly, we will recover this fea- 
ture in the turbulent (convective) dynamo described in 
section pi (see Fig. 13 'b).) 

This aTt dynamo provides a simple example of how our 
diagnostic may be interpreted in the context of stellar 
dynamo. Based on how our diagnostic highlights the 
saturating properties of the solar differential rotation in 
an a — ft case, we now apply it to a turbulent dynamo 
triggered in a stellar convection zone that also exhibit a 
solar-like differential rotation profile. 

4. NONLINEAR CONVECTIVE DYNAMO 

We use the general method described in Section [2] and 
validated in Section [3] to study dynamo action in a global 



(spherical) nonlinear convection zone. Contrary to Sec 
tion|3j we now solve the full set of MHD equations and do 
not introduce an y a effect. We model a turbulent solar 
convection zone dBrun et al. 2004 Jouve & Brun 2009 



Pinto fc Brun|2012[) that develops a solar-like different 1a 



rotation profile (.b ig. [4^a)) , with fast equator and slow 
poles. We display the convective patterns we obtain in 
Fig. gb). We recover the well-known 'banana'-shaped 
cells at the equator, and more patchy patterns at higher 
latitudes. Our choice of parameters yields a mildly tur- 
bulent state (based on the maximum amplitude of the 
velocity, the Reynolds number in the middle of the con- 
vection zone is of the order of 800). 

We display in Fig. [4jc) the kinetic energy spectra 
in the rotating frame at the center of the convection 
zone as a function of the shell L. We separate the 
axisymmetric component (the plain blue line) from the 
non-axisymmetric component (the dashed red line), 
and the dotted black line is the total spectrum. Notice 
that two peaks at L = 3, 5 dominate the kinetic energy 
spectrum. They represent the differential rotation of the 
azim uthal component of the toroidal velocity (see Sect. 
3.1). 



We initialize a peaked (/, m) = (9, 5) non-axisymmetric 
magnetic field (Fig. [5|a)) throughout the convection 
zone by setting: 



B = 



10B R% 



o5 
S 9 . 



(30) 



We set Bq = 100 G so the initial magnetic energy con- 
tained in the L = 9 shell is comparable to the kinetic 
energy at that scale (see Fig. [5|b)). 

The magnetic Prandtl number throughout the convec- 
tion zone is set to P m = ^ = 4, leading to a magnetic 
Reynolds number of the order of 3200 at mid- convection 
zone (based on the maximum amplitude of the velocity). 
Such a set of parameters triggers a dynamo inst ability 
and the growth of magnetic energy (see Fig. 6(c) in the 
following) . 

The initialization we chose allows us to directly see 
how a significant amount of energy can be transfered to 
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5.e + 02 nHz 



Fig. 4. — Progenitor hydrodynamical state, (a) Mean differential rotation profile averaged over 3 months, (b) Time-dependent convective 
patterns at the top of the convection zone with dark tones representing downflows. (c) Kinetic energy spectra in the middle of the convection 
zone. The axisymmetric spectra are in plain blue, the non-axisymmetric spectra in dashed red and the total spectra in dotted black. 




L shells 

Fig. 5. — (a) Initial magnetic energy (the seed poloidal streamfunction is a I = 9, m = 5 spherical harmonic) in cgs, in the middle of the 
convection zone, (b) Magnetic energy spectra one time step after the introduction of the peaked magnetic field in our turbulent convection 
zone. The axisymmetric spectra are in plain blue, the non-axisymmetric spectra in dashed red and the total spectra in dotted black, (c) 
Saturated magnetic energy in cgs, 600 days after the introduction of the seed magnetic field. The color scale is logarithmic. 



large scales. We also did the same numerical experiment 
varying the initial conditions. By initializing roughly the 
same amount of energy distributed over the whole scales, 
we obtained the same statistical saturated state. Hence, 
this proves that in this case, the initial scale is forgotten 
when the dynamo saturates. 

The complex interactions between the convective mo- 
tions and the initially peaked magnetic field lead to the 
construction of the magnetic energy spectrum. The sat- 
urated magnetic energy after 600 days of evolution is 
displayed on Fig. [5^c) in physical space, in the middle 
of the convection zone. In the remainder of this sec- 
tion, we characterize how such a state is obtained, and 
maintained. We distinguish two regimes: the develop- 
ment of the spectrum shape (the kinematic regime, Sect. 
4.1), and its sa tura tion and sustainment (the non-linear 
regime, Sects. 4.2 and |4.3| ). The results of this section 
will be summarized in Fig. 14 We recall here that no 



a effect has been added to tne induction equation ([5| 
dynamo action is n aturally achieved since convection is 
3D and R m > R™ {t ( |Brun et al.||2004 ). 



4.1. Creation of magnetic energy spectrum: 
kinematic phase 



We plot in Fig. 6(a) the evolution of the non- 
axisymmetric magnetic energy spectrum. The initial 
spectrum is plotted in blue, and the s atura ted spectrum 
in red. In addition, we display in Fig. 6(c)| the evolution 
of the magnetic energy for 6 different L shells. The total 
energy evolution is also shown (plain thick line). The 



initial evolution (t < 100 days) is shown in logarithmic 
scale The saturation of magnetic energy is reached at 
t ~ 300 days. 

We also ran another numerical experiment where we 
artificially suppressed the Lorentz force and the ohmic 
heating in the momentum and energy equations (i.e., ef- 
fectively running a kinematic dynamo). On average, the 
relative difference with the fully non-linear case starts 
being significantly different (departure of order one) 
roughly 10 days after the introduction of the magnetic 
field (the exact length of the kinematic phase depends on 
the scale considered) . We detail hereafter how the non- 
axisymmetric (Sect. 4.1.1) and the axisymmetric (Sect. 
4.1.2) spectra are created during these first days, which 



we will refer to as the kinematic phase. 

4.1.1. Creation of the non-axisymmetric spectrum 

We obse rve at first that all the L shells gain energy 
(Fig. 6(a)), excepts the L = 9 shell which looses energy 
because it is redistributed thro ughou t the whole domain 
by the convective flows (Fig. 6(c)). It stops decaying 
at t ~ 17 days. We identify lour regions in the non 
axisymmetric spectrum that exhibit different behaviors. 
We define the large-scale zone (I) by 1 < L < 4, the 
neighborhood zone (II) by 4 < L < 13, and the plateau 
zone (III) by 13 < L < 60. The small-scale zone (IV) 
(L > 60) starts at the highest diffusive scale, which is 
the highest viscous scale I ~ 60 based on the first scale 
at which the local Reynolds number is lower than 1. It 
also includes the magnetic dissipative scales (L > 120). 
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Fig. 6. — (a) Evolution of the non-axisymmetric magnetic energy 
spectra. The initial spectra is blue, the saturated spectrum red. 
Two slopes (L~ 6 and L~ 3 ) are given as references for the small 
scales, (b) Contributions to the non-axisymmetric magnetic energy 
evolution in the initial phase (dotted line in panel (a)). The total 
is in plain black, V in dotted black, V in blue and T in red. (c) 
Evolution of non-axisymmetric energy of shells 1, 3, 7, 9, 11 and 
42. Notice the fast early evolution of E% and E*-^ due to shearing 
of Bg by the differential rotation. The thick plain line is the total 
non-axisymmetric energy. The abscissa is in log scale between for 
t < 100 days, and in linear scale for t > 100 days. Even though 
the initial growth is exponential, we chose to represent it in log-log 
scale to make it appear clearly in the evolution plot. 



The four zones a re sep arated by the three dotted verti- 
cal lines in Fig. |6(a)| In order to understand how the 
spectrum is built, we di splay t he co ntributi ons fr om the 
different terms of Eqs. ( |A25| >-( |A28| > in Fig. |6(b)| Those 
contributions are taken shortly alter the introduction of 
the magnetic field. They corresp ond t o the spectrum 
plotted with a dotted line in Fig. |6(a)| We recall that 
we fully calculate all the coupling terms up to L = 70. 

The energy transfers around the L = 9 shell (the neigh- 
borhood zone II) are dominated by U° — B* interactions 
from both the production V and the flux T terms. We 
recall here that both V and T represent generic transfer 
functions, that can either be positive or negative. Dissi- 
pation is negligible in zone (II), even for the L = 9 shell 
that initially contains the energy. The L = 9 energy then 
decreases through the interaction of Bg and the differ- 
ential rotation U3 that shears the magnetic field (see 
Sect. 3.1). The energy is preferentially redistributed to 
E\ x and Ej. For those two shells, the production V and 
flux T terms co ntrib ute positively to the creation of the 
spectrum (Fig. |6(b)[ ). We display in Fig. [7] the detailed 
contribution of P and T to Ej. We only display con- 
tributions from B* because the axisymmetric magnetic 
energy is very small initially. The U* — B* interactions 
are displayed in panel (a), and the U° — B* interactions 
in panel (b). We sum over the velocity shells to plot the 
production term against B* in panel (c). We observe 
that the summed contribution is dominated by Bg — U3 
interactions, as expected. Also, we observe that energy 
is directly transfered from Eg to Ej , such that the L — 8 
shell is not involved in the transfer. This is true for all 
the shells in zone (II) and implies that the transfer of 
energy is non-local, even for shells close to the initial 
energetic shell. 

Due to the triangular selection rule, the U3 — Bg in- 
teraction can only act in zone (II). Indeed, L must be 
strictly greater than 12 in zones (III-IV) and strictly 
lower than 6 in zone (I). U3 and Bg initially domi- 
nate respectively the kinetic and magnetic energy spec- 
tra. Their interaction was consequently dominant in zone 
(II), and we expect a different kind of spectral transfers 
in the other zones. This zone exists because of our choice 
of initial condition. The very early evolution would have 
been changed if we had chosen a different initial shell. 
Though, as stated before, this initial sc ale is forgotten 
when the saturated state is reached (Fig. 6(a)). 

The dynamics of zones (I), (III) and (IV) are domi- 
nated by two effects which competes initially: a direct 
non-local Bg — U* transfer of energy, and an effective 
shearing of neighbor shells by the large scale differential 
rotation (B* — U3 interactions). These two effect are 
exemplified in Fig. [8] for L = 42 (zone III). 

In the case of the Targe-scale zone (I) (not shown here), 
the evolution is dominated by both V and T . The in- 
teractions between Bg and U* alternate signs depending 
on the U* shell considered. We also stress that the in- 
teractions involving other B* shells are not negligible. 
The differential rotation action is completely negligible 
compared to U* — B* interactions in zone (I). 

In the case of the plateau zone (III), almo st a fl at pro- 
file in the log-log plot is observed in Fig. |6(a)| (hence 
its name). This plateau is characteristic or connective 
flows that usually exhibit a broad spectrum between the 
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Fig. 7. — Production V and flux T contributions to the non-axisymmetric shell L = 7. The 2D color maps are the B* — U* and B* — U° 
transfer functions, the ID plot is the sum of the transfer functions over the U shells. The horizontal dotted line labels the B7 shell. The 
oblique black lines represent the boundaries of the triangular selection rule. 
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Fig. 8. — Production V contribution to the non-axisymmetric shell L = 42 



injection and inertial ranges (Fig. (4^c)). The evolution 
of the spectrum is dominated only by the Vl (Li, L2) 
contributions (F is negligible), and in particular by the 
coup ling between (non-axisymmetric) U* and B* (Fig. 
6(b)). Hence, it is a non-local transfer of magnetic en- 
ergy that creates the spectrum. All the shells in zone 
(III) receive energy mainly through this non-local mech- 
anism. As a result, the energy transfer is very sensitive 
to the kinetic energy contained in the U£ shells involved 
in the coupling. This explains why the magnetic energy 
spectrum reflects the kinetic energy spectrum in this re- 
gion. 

Although the U* — B* interactions dominate (Fig. 
|8j, we stress that the U° — B* interactions exhibit a 
airect cascade pattern. E\ 2 receives energy from E% Q 
through U3 — B4Q interactions, and gives energy to E\± 



through U3 — B^ interactions (see panel (b) in Fig. 
[8|. Even if the triadic interaction involves the large 
scale velocity U3, we nonetheless refer this effect as a 
cascade. The velocity field only acts here as a mediator, 
and the scales of magnetic field involved in the magnetic 
energy transfer are at the same scale. It is consequently 
a cascade when considering the scales of magnetic field. 

The energy transfers in zone (IV) (not shown here) 
are very similar to zone (III). A noticeable difference is 
that the cascade of energy triggered by the shear of the 
differential rotation U3 is much less efficient since the 
smallest scales hardly feel the large scale rotation pro- 
file. Finally, ohmic diffusion acts in the whole zone (IV) 
and tends to dissipate energy. It has a sufficiently lower 
amplitude than the non-local transfers so that it does not 
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dictate the spectrum shape initially. It wil l nev ertheless 
contribute to the saturation process (Sect. |4.2[). 
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9. — Same as Fig. [6] for the axisymmetric part of the spec- 
In panel (a), the initial spectrum is zero. 



We now characterize the cre ation of the axisymmetric 
spectrum. We display in Fig. |9(a)| the evolution of the 
axisymmetric component of the magnetic energy. We re- 
call that since we initialize the dynamo with a purely 
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Fig. 10. — Production V for the axisymmetric shell L = 1 during 
the initial state. Both U° — B° and U* — B* couplings are shown. 



non-axisymmetric field, the initial axisymmetric spec- 
trum is null. After one time-step, the axisymmetric mag- 
netic energy is orders of magnitude lower than the non- 
axisymmetric spectrum (Fig. j5^b) The global shape of 
the axisymmetric spectrum is created very rapidly, all 
the shells gain energy at about the same rate until they 
saturate. The initial exponential growth rate is the same 
for both the axisymmetric and non-axisymmetric spec- 
trum is approximately 0.6 days -1 (which corresponds 
to a time-scale approximately 17 times lower than the 
convecti ye tu rn-over time). This can also be observed 
on Fig. 9(c)| where we plot the evolution of few shells 
against time. They all gain energy at about the same 
rate initially, and then slowly tend to a saturated state. 
The axisymmetric shells considered have comparable 
energy since the spectrum is essentially flat at scales 
L < 30 (Fig. [9(a)] ), which was not the case for the non- 
axisy mmetric spectrum (Fig. |6(c)[ ). We observe in Fig. 
9(b)| that the flux term T plays a major role between 
L = 2 and L = 13. This means that the creation of the 
spectrum is dominated by the radial interactions at those 
scales. The two flux curves exhibit a sawtooth pattern 
that is again reminiscent from th e diffe rential rotation 
energy shells (see Sects. |3.1| and 



4.1.1). At higher L, 



the evolution of the spectrum is the result of a complex 
interplay between the production and flux terms. 

More interesting, the dipole {L = 1) evolution is 
dominated by the production term through the interac- 
tion between the non axisymmetric magnetic field and 
velocity field. We display the detailed transfers maps 



for this scale on Fig. 10 We observe that the large scale 
magnetic field is mainly created by the interplay between 
By and Ug_ 8 . The transfers involving Bg (where the 
energy is originally mainly contained) act negatively and 
do not dominate the transfer of magnetic energy. This 
is consistent with the fact that the whole axisymmetric 
spectrum shape is rapidly created and only gains energy 
globally afterwards. It does not depend on the scale at 
which we initially put the non-axisymmetric magnetic 
energy. Since the energy is not transfered directly 
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from the initial reservoir of energy .Eg, we already 
see preferred transfers towards the large scale dipole 
involving B£, which is one of the highest energy scale 
of the non-axisymmetric spectrum at this time. This 
effect shal l be confirmed during the saturation phase 
(Sect. 4.3). The creation of the axisymmetric magnetic 
energy spectrum seems to depend essentially on the 
initial hydrodynamic convective spectrum (as expected 
in such kinematic phase). 

4.2. Non-linear saturation of the smallest scales 
Following Sect. |4.1[ we now detail the saturation 



the horizontal dashed line [L = 152), and by negative 
contribution (blue) above. This is confirmed by the 
plot in panel b where the transfers have been summed 
over the velocity shells. This cascade is of different 
kind than the one observed in zone (III) (Fig. pTj). 



and sustainment of the magnetic energy spectrum at 
small scales. By 500 days the axisymm etric and non- 
axisymmetric spectra are saturated (Figs. |9(c)| and|i 

y to never be 



- |6p 
De null at 



The flux T contribution is likely 
the largest scales since it represents the flux of magnetic 
energy through the horizontal surface at the middle of 
the turbulent convection zone. In order to saturate the 
magnetic energy (i.e., to get dE™ ag /dt = 0), V and/or V 
have to compensate T . In the first three zones, diffusion 
is negligible. Hence, V natur ally tends to cancel T out 
in those zones (see Sect. |3.2| for a simple version of this 
cancellation effect). The cancellation effect is such that 
.F(U*,B*) tends to cancel P(U*,B*) out. This is also 
the case for J*(U°,B*), J*(U*,B°), and J*(U°,B°). 

In spite of the cancellation of the different contribu- 
tions, characteristic patterns can still be identifi ed. The 
more distinctive pattern we identified in Sect. |4.1| was 
the direct cascade of magnetic energy in zone (III) . It 
turns out that we still observe it and that it slightly dom- 
inates the transfer terms during the saturation phase. 
We display on Fig. [IT] the production contribution to the 
non-axisymmetric magnetic energy evolution 600 days 
after the magnetic field introduction. We recover the 
direct cascade of energy in the production contribution, 
that was already present on Fig. [8] This direct cascade 
of energy is associated with an inverse cascade of energy 
carried by the flux contribution, which opposes the 
production term during the saturation phase. Both 
cascades are of the same order of magnitude and tend 
to cancel each other out. They are associated with the 
axisymmetric component U3 (the differential rotation), 
and the non-axisymmetric components of B. The 
contributions of non-axisymmetric components of U 
involve more shells, but their net effect is a bit lower 
than the shear from differential rotation (panel (a) on 
Fig 11). On this panel, no particular global pattern can 
be ioentified. 

The transfers of magnetic energy appear to be very 
interesting in zone (IV) where diffusion acts signifi- 
cantly. In order to saturate, V and T have to combine 
to cancel V. For the non-axisymmetric spectrum, it is 
the production term that dominates over the flux term 
to compensate diffusion. In addition, the production 
term in zone (IV) exhibits a very particular generalized 
cascade shape. This cascade could not be identified 
during the early evolution for it was dominated by the 
non-local transfer from Bg. We display on Fig. 12 (panel 
a) the U* — B* production map towards Ef^. The 
other interactions are negligible. We observe that the 
map is dominated by positive contribution (red) under 



Here, no clear velocity shell dominates the transfer map 
(panel a on Fig. [l2| ). It is a generalized cascade that 
results from the coupling between many magnetic shells 
(around L = 152) and all the largest velocity scales. 
Hence, the velocity scales involved in the cascade are 
not local compared to the magnetic field scale considered. 

Trying to simplify the complex 2D transfer maps, 
one may isolate the main contributing couplings to the 
different evolution terms. Doing so at all times for the 
non-axisymmetric spectrum at small scales, we find that 
the percentage of couplings that account for 90% of the 
contributing terms typically varies from nearly 1 to 70% 
of the calculated couples. As a result, we demonstrate 
here that the complex dynamo process occurring in 
a 3D turbulent convection zone involves many modes 
that interact though non-trivial triadic interactions. 
Then, the dynamics of the smallest scales can hardly be 
reduced to the evolution of a small set of modes. 

Finall y, th e analysis of the axisymmetric aQ dynamo 
in Sect. |3.2| shed light on the importance of the families 
of symmetry (with respect to the equator) of the fields. 
The instantaneous convective motions do not exhibit any 
particular symmetry at any scale and the kinetic energy 
spectrum is a mixture of both primary and secondary 
velocities. The differential rotation is the only velocity 
featur e th at has a clear symmetry (secondary family, see 
Sect. 3.1) and that has a large influence on the mag- 
netic energy spectrum. It is involved in the magnetic 
energy cascade in zone (III), and shears both primary 
and secondary magnetic fields to cascade primary and 
secondary magnetic energy. Thus, it does not select a 
particular symmetry. Indeed, the ratio of primary (an- 
tisymmetric) to secondary (symmetric) magnetic energy 
varies with time for all shells and does not settle even 
during the saturation phase. The presence of complex 
flows, often breaking the equatorial symmetry, yields a 
str ong coupling o f both dynamo families (as in the Sun, 
see |DeRosa et al.| ( |2012[ )), contrary to simpler mean field 
dynamo models (see Sect. |3|. 

4.3. Sustainment of the mean large scale magnetic field 

Given their key r ole in setting the o verall magnetic 
polarity in the Sun ( DeRosa et al.|[2012[ ), we now detail 
the main contributions to the saturation and sustainment 
of the large scale axisymmetric dipole (Z = 1, m = 0) and 
quadrupole (Z = 2, m = 0) fields. 
At the late phase of the simulation the lar ge-sca le ax- 



isymmetric spectrum is fully saturated (Fig. 9(c)). The 
saturation is obtained thanks to the compensation of 
the production and flux terms, similarly to the satura- 
tion of the mid-scales (see previous section). The large- 
scale dipole (I = l,m = 0) saturation process differs 
significantly from its creation. We display in Fig. [13] 
the production maps for the axisymmetric dipole and 
quadrupole averaged over 150 days during the saturated 
state. The transfers maps of T (not shown here) are 
exactly opposite to the maps (a) and (c) for V. We 
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Fig. 11. — Production V contribution to the non-axisymmetric shell L = 42 during the saturation phase. Interactions between B° and 
U* are negligible. 
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Fig. 12. — Production V contribution to the non-axisymmetric 
shell L = 152 during the saturation phase. The red line is the 
smoothed contribution, that exhibits the characteristic S-shape of 
the generalized cascade. 



see that both the axisymmetric and non-axisymmetric 
fields significantly contribute to the saturation and sus- 
tainment of the large scale dipole. In particular, two 
main contributors emerge. First (panel a), the coupling 



of the differential rotation U? with the large scale B4 



field dominates the axisymmetric contributions. This ef- 
fect is more likely to represent the shearing of the large 
scale poloidal multipole by the large-scale toroidal 
differential rotation. 

Second, the non-axisymmetric contributions (panel c) 
are at least equally important for the saturation of the 
dipole. In particular, the interaction U23 — B23 domi- 
nates the non-axisymmetric contributions. Thus, is it a 
non-local interaction that saturates the large-scale mag 
netic dipole. Furthermore, B23 is one of the most ener 



getic shells of the magnetic energy spectrum (Fig. 6(a) ). 
This directly points out the importance of the mid-scale 
part of both the kinetic and magnetic energy spectra for 
the saturation level of the large-scale magnetic dipole. 

We can remark here that the major contributions of 
V for the saturation of the dipole are all positive. They 
are balanced by negative contributions from T . Conse- 
quently, if the differential rotation was more efficient, or 
if the U23 — B23 interaction possessed more energy, the 
saturation level of the large scale dipole would be much 
higher. 

Since our magnetic Prandtl number is 4, the peak of 
the kinetic and magnetic energy spectra are likely to 
be shifted. At saturation, the couplings are nonetheless 
dominated by the peak of the magnetic energy spectrum 
that occurs at smaller scale than the peak of the kinetic 
energy spectrum. Changing the magnetic Prandtl num- 
ber will cause the separation of the peaks to change. If 
the peaks separate more, our results suggest that the 
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Fig. 13. — Saturation of the large scale axisymmetric dipole (pan- 
els a and c, production V for E®) and quadrupole (panels b and 
d, production V for E®). The interaction maps are time averaged 
over a period of ~ 150 days. Both U° — B° and U* — B* couplings 
are shown. 



saturating interaction will involve smaller scales velocity 
and magnetic fields. The velocity field involved is likely 
to be less energetic, which could trigger a smaller saturat- 
ing interaction, and in turn a lower energy state for the 
large scale dipole. If the peaks are closer (or eventually 
switch), the picture becomes more complicated and we 
cannot predict if the saturating interaction will remain 
fixed by the peak of the magnetic energy spectrum. The 
exploration of this parameter space is left for future work. 

The large scale quadrupole also saturates thanks 
to both the axisymmetric and non-axisymmetric fields 
(panels b and d) . The axisymmetric contributions (panel 
b) are very similar to the dipole case and are domi- 
nated by the differential rotation. The differential ro- 
tation shears both and B5, which is opposed by the 
flux term to saturate the quadrupole. Again, this effect 
accounts for the saturation of the toroidal quadrupolar 
field. Hence, the saturated level of the poloidal dipolar 
field (panel c) plays a major role for the saturation of the 
toroidal quadrupolar field. 

The poloidal quadrupolar field is then saturated 
through the non-axisymmetric interactions (panel d). 
The contribution are again very non-local, though in this 
case no particular scale dominates the saturation pro- 
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cess. Hence, we may expect that the saturation process 
of the axisymmetric quadrupole will have a very differ- 
ent dependency on the magnetic Prandtl number than 
the axisymmetric dipole. 

5. CONCLUSIONS AND PERSPECTIVES 

In this paper we developed and validated a new spec- 
tral analysis method suited for spherical objects. Using 
two vectorial spherical harmonics basis, we were able to 
calculate transfer functions of magnetic energy in spec- 
tral space. We can calculate the coupling coefficients up 
to I ~ 500. For the first time in such studies, the com- 
plete 2D transfers maps have been calculated to charac- 
terize the full triadic interactions. 

After a quick numerical validation, we first applied 
our method to a simplified dynamo case. Such 
axisymmetric m odels are very wel l known to trigger 
cyclic dynamos ( jCharbonneau 2010| ) with our choice of 
a symmetric (with respect to equator) velocity field and 
an antisymmetric a effect. The clear separation between 
the dipolar and quadrupolar families was illustrated 
thanks to our new diagnostic. The production (i.e., on 
a spherical surface) and a flux (i.e., through a spherical 
surface) contributions were shown to quasi-cancel each 
other out for all shells. 

Our method was then successfully applied to a 3D 
turbulent convective dynamo case. We initialize a highly 
non-axisymmetric magnetic field and let the dynamo 
develop a turbulent spectrum of magnetic energy. We 
distinguished the kinematic phase with exponential 
growth of the magnetic energy spectrum, and the non- 
linear ly saturated phase. The first phase is dominated 
by both a non-local transfer of energy from the initial 
scale of magnetic energy, coupled with the convective 
scales, towards all the other magnetic scales, and the 
shearing by the large scale differential rotation. A large 
part of the magnetic energy spectrum is then dictated by 
the kinetic energy spectrum developed by the convection. 

The saturation phase is more subtle and greatly de- 
pends on the considered scale in the spectrum. The sat- 
urating interactions for the different spectral scales are 
illustrated in Fig. [l4j and summarized hereafter. 

Our new methodallowed us to distinguish two clear 
cascades of magnetic energy at the smallest scales of our 
simulation, for 13 < L < 60 and L > 60 (highest L's). In 
the former case, the differential rotation profile mediates 
the cascade by shearing the magnetic field. It results in 
an efficient cascade of magnetic energy. 

The latter cascade is also direct and involves all the 
highest velocity scales (the large scale differential rota- 
tion does not dominate in this case). It is a generalized 
cascade over a large range of magnetic scales. The 
velocity scales involved in the cascade are not local with 
respect to the magnetic scales. As consequence, we 
cannot predict if this generalized cascade would hold 
at the lowest scales in the case of a real convective 
dynamo where scale separation is much higher. Besides, 
the saturation also involve non-local coupling that can 
eventually be of the order of the sheared cascade for 
the intermediate scales. We proved in that case that 
the transfers cannot be reduced to a limited set of modes. 



The saturation of the large scale axisymmetric dipole 
and quadrupole appear to be radically different than 
the small-scale saturation. The toroidal components are 
mainly saturated by the balance of the shearing effect 
of the differential rotation on the large scale poloidal 
fields, and the flux transfer through the spherical shell 
due to the effect of the differential rotation. The 
poloidal components are mainly saturated by non-local 
non-axisymmetric interactions. The dipole is saturated 
by the scale of maximum (highest) magnetic energy, and 
the quadrupole saturation is not dominated by any par- 
ticular scale. These two observations point to the two 
main dependencies of the saturating interactions for the 
large-scale fields. First, the rotation rate of the star 
(which is linked to the saturating interaction through the 
differential rotation) ca n determine the ab ility of the dy- 
namo to build wreaths ( |Brown et al.|2010|), and/or to be 
in a strong or weak regime ([C hristensen & Aubert 2006 



Featherstone et al.|2009||Simitev ^"Bu sse 2009). Second 



the magnetic Prandtl number F m = vj r\ determines the 
postion of the peak of magnetic energ y and then affects 
the s aturating interaction (e.g., see |Schekochihin et al. 
2004). We will explore in detail how the saturating in- 
teractions depend on those two effects in future work. 

Finally, it is worth comparing theses resul ts with 
previous related work of Livermore et al. (2010). When 
using forced helical flows and allowing the dynamo field 
to back-react on the flow, they found that the saturation 
of the large scale poloidal dipole was dominated by 
non-local interaction with a particular magnetic scale 
(I = 10, m = 10 in their case). They showed that 
magnetic energy was transfered to this scale by the 
large-scale toroidal magnetic field. In our case, the 
large-scale dipole is also saturated due to non-local 
interactions. Though, the flow we consider is signif- 
icantly different because (i) it is obtained from the 
convective instability and (ii) its spectrum is dominated 
by the large-scale differential rotation that develops 
self-consistently. Hence, the dynamo process is different 
and we find that the toroidal large-scale field is saturated 
by the effect of the large scale differential rotation, and 
the large-scale poloidal dipole by the non-local transfer 
of energy from the magnetic scale of maximal energy. 

Our results also suggest that no significant large-scale 
magnetic field is growing over dissipative time scales in 
our simulation (the ohmic dissipation time scale for the 
axisymmetric dipole is typically of the order of t v ~ 700 
days in the simulation). Again, the fast saturation of 
the dynamo (less than 300 days) may not hold for lower 
magnetic Prandtl number dynamos. 

We developed a diagnostic on the magnetic energy 
that is an invariant of ideal MHD. In the case of 
non-ideal MHD, the existence of the selective decay 
dTaylor|[l 974; Mattha eus" fc Montgomery ||1980| |Mininni| 
|&; Mo ntgomery 2006) introduces a decoupling between, 
e.g., the evolution time scales of the magnetic (or total) 
energy (fast) and the magnetic helicity (slow). As 
mentioned before, we were interested, in this work, in 
fast phenomena compared to the ohmic diffusion time. 
For such processes, the ideal invariants of MHD are 
still the appropriate quantities to interpret the scales 
interactions. The dynamo saturation is necessarily 
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Fig. 14. — Saturating interactions of the non-linear convective dynamo summarizing the results of section^] KE and ME stand for 
Kinetic and Magnetic Energy and DR for Differential Rotation. 



achieved through a modification of the kinetic energy 
spectrum. As mentioned before, the case we studied in 
this paper is in the weak branch of the dynamo (i.e., the 
large scale poloidal magnetic fi eld does not dominate 



the magnetic energy spectrum |Christensen fc Aubert 
2006; |Simitev fc Busse||2009| |Gastine et al.||2012| ). Tne 



detailed modification ol the kinetic energy spectrum 
implied by the saturation of the dynamo process will 
be discussed for both the strong and weak branches 
in a future work. Although the magnetic energy is 
the relevant quantity to characterize nearly kinematic 
dynamos (where the Lorentz force plays little role), a 
diagnostic on kinetic energy would be highly valuable 
for non-linear ly saturated dynamos. On top of that, 
the detailed spectral transfers of magnetic helicity are 
also mandatory to fully address the complexity of the 
dynamo process. Evolution equations of kinetic energy, 
magnetic helicity and cross helicity in the framework 
introduced in this paper are under development and will 
be published in a forthcoming paper. 



Finally, those diagnostics may also prove very useful 
for non-dynamo transfers related MHD phenomena. For 
example, spectral analysis applied to the relaxation and 
the stability of low-1 fossil field (see |Braithwaite fc Nor d- 
hmd|2006l ^nm|f2007l [Zahn et al.||2007[ |Duez fc CTatT 
2010; IDuez et aT|20To|) will be studied in a future publi- 



cation. 
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APPENDIX 

A. DEFINITION AND PROPERTIES OF VECTORIAL SPHERICAL HARMONICS 
A.l. Classical vectorial spherical harmonics basis 
A. 1.1. Definitions 
We define from |Rieutordl fl 198 7Q ; |Mathis fc Zarm| ( |2QQ5[ ): 

S?(0,<p) = V ± Yr =d e Y l m e e 
T?(0,<p) = V ± x Rr = s^lfe, - W% 

where (e r ,eg,e v ) defines the spherical basis and F ; m are the Laplace spherical harmonics defined by 



Y, r > 



,<p) = ("!)■ 



1(21 + 1) (l-H)!„ H 



4tt (l+\m\)\ 1 



Pr l (cos0)e mv 



where P™ are the associated Legendre polynomials. The basis ([AT]) have the following properties : 

f ^-(R^r dn=s lul2 s mum2 , 

Js 

[ s^.(S™ 2 ) cc dQ=f T^.(T™ 2 ) CC dn = h(h + l)5 luh 5 mum „ 
Js Js 



(Al) 



(A2) 

(A3) 
(A4) 



where S is a spherical surface, = smOdOdcp the solid angle, cc means complex conjugate and S is the Kronecker 
symbol. We also have: 



(spr = (-irs,- 



(A5) 



and all the other scalar cross products are 0. We remind the reader that the poloidal fields are described by their 
projection on (R] 71 , S] 71 ), and the toroidal fields by their projection on T] 71 . 
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A. 1.2. Scalar fields 
Defining i/>(r,0,<p) = E^i EL=-z {^m(O^ m (0> we S et: 

oo I 
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^=E E { 

1 = 1 m=-l k 
oo I 

v-vv>=e E A ^ 



(A6) 
(A7) 



1 = 1 



where A t = d 2 rr + \d r 



A. 1.3. Vectorial fields 

For a vector X(r,0,<p) = Ej=i EL=-j {-0^ + W + (r)T™ }, we obtain: 

OO £ 

v-x=EE 



1=1 rn = — l 
oo Z 



^a r (r 2 ^)-Z(Z + l)^ 



VxX = ^ E 

£=1 m = — 

oo Z 

v 2 x=E E 



Z(Z + 1)- 



r? 



d r (r<) 



£=1 m = — l 



AiB l m +2 



^A l m --,(A l m -l(l + l)B l m ) 



R7 



.4' 1 



(A8) 
(A9) 



A 1 



Sr + [A ; C^]T m 



(A10) 



A. 1.4. Recurrence relations 



In addition to the expression of the di fferent operators, we also give here two useful coupling relations between 
spherical harmonics. First, according to Varshalovich et al. ( |1988| ), the coupling between cos# and the spherical 
harmonic Y™ is given by 



cos ™ _ 4 (l-m+l)(l + m+l) 
C0S6Yl "Y (2Z + l)(2/ + 3) Yl+1 

Then, one easily deduces the following properties: 



' (I — m) (I + m) 
(21 - 1) (21 + 1) 



r i-i - 



(All) 



cos<9S™ = 



I (l-m+l)(l + m+l) 



Z(Z + 1) 1 ' 



(21 + 1) (2/ + 3) 



'i+i 



Z_+l (l-m)(l + m) 

~y (2i- 1) (2i + i) 



Sm 
i-1 



cos6»Tj™ = 



Z /(Z-m+l)(Z + m+l) 



z + i V 

_jrn_ 
I (* + !)' 



(2Z + 1) (2Z + 3) 



A H-1 



Z + 1 (l-m)(l + m) 
TV (2Z — 1) (2Z + 1) 



rr-im 
A J-1 



(A12) 



(A13) 



A. 2. alternative vectorial basis 
A.2.1. Definitions 

The vectorial spherical harmonics basis defined in appendix |A.1| is very efficient to calculate scalar products or linear 
differential operator on vect ors. Nevertheless, it is quit e hard to use it to express vectorial products. Instead we define 
the following basis (e.g., see Varshalovich et al. ( 1988[ )): 



Y l J +v {0,<p) = E {(-l)'- m v / 2Z + l( 
n=-i ^ ^ 



I l + v 1 
m i± — rn — /i I 1 1 +" C/x ( ' 



(A14) 
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where (...) is the 3-j Wigner coefficient linked to Clebsch-Gordan coefficients, and the vectors e M are 



eo = e z 



where (e x ,e y ,e z ) defines the cartesian basis. Note that the equivalent of the conjugation rule (A5) is then 



\ X l,l+v) — X l,l+v ' 



(A15) 



(A16) 



Again, we recall that the poloidal fields are described by their projection on (Y[ r J +1 , Y[ r J_ 1 ) (y £ { — 1; 1}), and the 
toroidal fields are described by their projection on Y™ (y = 0). 

A. 2. 2. Vectorial product 
We decompose a vector X on this basis in the following way: 

CO I 1 

x(r,e,<p) = Yl E E xM) Y i7 + v 

1=1 m = — l v— — \ 

Evaluating the vectorial product of two vectors X12 = Xi x X2, one gets: 



y m i2 

^12^12,^12+^12 



E\ A \ ^ Y 1711 Y 777,2 /t^ 12 ' 77712 '^ 12 

/ j / j ^l;h,h+vi 2\h,l2+V2' J li,mi,i>i,l2,r, 



h,l 2 =l m 1 = — l 1 vi,f2 
li2>\h-h\ rn 2 = -h 
li2<h+h m 1 +m 2 =mi2 



where 



l\ , m l i v \ ^2 i^l i v 2 ^ ' 



— (2Zi + l)(2Zi + 2 Vl + 1)(2Z 2 + 1)(2Z 2 + 2v 2 + 1)(2Z + 1)(2Z + 2^ + 1) 

Z7T 

il h I * 

h + v\h + v 2 l + v 
1 1 1 



h h I \ fh + v\ h + ^2 / + v 

mi rri2 —(mi + 777-2) / I 



with {• • • } being the 9-j Wigner coefficient. 



A. 3. Basis change relations 
For a vector X decomposed in the following manner: 

00 1 

x = E E M™ R r + & m s? + c L t t } 

1 = 1 rn = — l 

co I 1 

£=1 m=—l v— — l 

we have the two following relations to change from one basis to the other: 
A 1 - 



/2Z+T 

1 



1,1+1 



V27+T 

i y 77 



1 y m 



1 y m 



LZ-1 - 



S i 3 (^ + (Z + 1)B{ I 



L f — Al 4. /Hi 
1 V -^m ^ "-'m 



i+i 
2i + 



(A17) 



(A18) 



(A19) 



A. 4. Expression of the a effect 

The a effect introduces the spectral coupling of a scalar and a vector, which was not treated before. In the special 
case of an axisymmetric a and an axisymmetric vectorial field B^e^ (which is the case in this paper, see Sect, pi), we 
write the a coefficient 

+00 

a(r,e) = Y,ao(r)Y l °, 

1=0 
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and we rewrite the magnetic field from (11) 



Introducing the coefficient 



— A 1 



1=0 



1(1 + 1)Z 2 (Z 2 + 1)(2Z + l)(2Zi + 1)(2Z 2 + 1) 



47T 



we can write the a effect such as 



h k I 




h k I 
1 -1 



(A20) 



aB^e^ 



+oo 

E 

1=0 



Finally, one gets 



h,h=o 
l>\h-h\ 

v l<h+l 2 



+ 00 W 

E| T ?- 



(A21) 



Z=0 



V x Ke,) = ^ + \d r @) S?} . 

1=0 ^ ' 

A. 5. Couplings in the magnetic energy equation 



(A22) 



The detailed expressions of the different terms of the magnetic energy equation (17)-(19) are given here. We write 
the magnetic field B and the current J 



+oo I 1 

B = E E E B^(r)Y^, 

1 = 1 m=—l v= — \ 
+oo I 1 

J = E E E ^+» y m+,- 

1 = 1 m= — l v= — \ 

In this basis, the vectorial product may be evaluated thanks to a coupling coefficient J[ 

(A18). The transformation rules from one basis to the other are given in appendix A. 3, and they allow us to easily 
evaluate the integrals (17)-(19). By separating the diffusive terms into two V\ and P 2 terms, we get that 



l,m,v 
lj ,m i ,vi ,l 2 ,m 2 ,V2 



(A23) 
(A24) 

given in equation 



P 1 (L,r)= ?7 ^(-irZ(Z + l) 

L 

cd r rj 



1(1 + 1) 



rA, 



E>T td m j-m y ' ' 1 
/ , n l,l+vi J l,l+is 2 J l,m,vi,l,- 



■m,V2 ' 



, (A25) 
(A26) 



where stands for a summation over all the spherical harmonics contained in the shell L (one element in an 
axisymmetric shell, and 21 elements in a non-axisymmetric shell). The production and flux terms then read 



T>L{L u L 2 ,r) = ±Y, 



(-l) m 1(1 + 1) 
r 2 ^/W+l 



^ (U Ll x B ia ) ™_ 1 +-^= (U Ll x B L2 ) - +1 



VI (U Ll x B L2 ) -VZTI(U L1 x B L2 )™ +1 



/ 2? + l 



(U Ll x B L2 )Jr 2 A, 



7- i (L 1 ,i 2) r) = -^aJr 2 ^ £ £ (U Ll xB i2 ) 



4 



0,0,1 



l,l-\-vi l,l+V2 *^l,m,vi ,1, — m,V2 



L V\=—\V2- 



(A27) 
(A28) 



The laplacian formula used for T>\ in the vectorial spherical ha rmon ics basis is given in equation (A10). In the 
production term we simply made use of the basis transformation (A19). Finally, the expressions f or X> 2 and the flux 
term need some intermediate steps to be properly explained. These details are given in Appendix |A.6| for the flux of 
magnetic energy, and the same procedure may be applied in the case of the second diffusive term. 
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A. 6. Simplification of the magnetic energy flux 
The flux of magnetic energy can be simplified, if one notes that it has the general form 

r. OO I 

F = / V • X dn and X = £ £ {^R^ +^1" + 07} . 
Then, one can easily deduce that 

F = [ V • X dtt = ^d r (r 2 A° ) . 
J s ^* 



Recalling from the system (A19) that R[j = — Y[j l5 and if one assumes that X = Xi x X2, one obtains, for an integral 



similar to the magnetic energy flux: 



JS I l=0m=-lv 1 ,v 2 ) 

A. 7. On primary and secondary families 

Previous studies on dynamos in stars shed light on the important distinction of primary (or dipolar, antisymmetric 
with respect to the equator) and secondary (or quadrupolar, symmetr ic with respect to the eq uator) families of magnetic 



field. For a vector X = YT=i TL=-i {^m^T + Bl m S T + C in T T}^ Roberts & Stix (1972) define the primary family 

as 

^ —^m ^ra+l 1 °rn ' u m m "~ ^ra+S ' °m ' u ra J -ra+2 ~T • • • 

and the secondary family as 

vs /im-pm j_ K?mom j_ /->ra+lrpra j_ yjm+2-nm j_ wm+2om j_ /->ra+3rpra j_ 

^ — ^m^m i" °m^m u m 1 m+1 rt m+2 ^ °m ^ra+2 ^ u ra - 1 m+3 • • • 

It can also be easily shown that in the ^YJJ] +I/ ^ basis, a primary field always satisfies I + m + ^ even, and a 

secondary field always satisfies I + m + ^ odd. 
Note that the vectorial product depends on the 3-j Wigner 

'i + v\ h + ^2 ' + ^ 


Recalling that mi+m2 = m, this 3-j Wigner is zero if /i+^i+/2+^2+^+^ = (^l+^i+miJ+^+^+^J + t^+^+^-^m 
is odd. Consequently, in order to have a non-zero 3-j Wigner, if U and B are from different families, their vectorial 
product is a secondary field; and if they are from the same family, their vectorial product is a primary field. If 
C = A x B, this means that 

^ff h C! and A* x W ^ CP ' (A29) 



where the superscripts p and s stand for primary and secondary. This was already acknowledged by [McFadde n et al. 
fl!99ip and [Gubbins fc Zhang| fl!993| ). 

B. NUMERICAL VALIDATION 

In order to validate the way we implemented in the ASH code the complex interactions between spherical harmonics, 
we compared an analytic calculation for a simple setup with numerical results. We summarize here those calculations. 

We start from a mixed m = 0, ±1 and I = 1 state for the magnetic field, and an (Z, m) = (2, ±1) state for the velocity 
field. We initialize the magnetic field in the following way: 

B (rs + i R! - Ir,-) - (s; + is; - is,-) + ^ (t; + ix; - i T -) . 

where Rq is the solar radius, is our inner boundary radius, R t is our outer boundary radius, and j3 = (R t — r) 2 (r — 
Rb) 2 - The velocity is initialized by: 

- ^ (Ri - RJ1) + m (s >- s;1) + *k {Tl - T; '» • 
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TABLE 1 

Analytical and numerical values of the vec tor ial product U x B (validation case, 

see Sect. |2.3|). 



SH Mode Analytical Expression 



Analytical Value 



Code Output 



(1.0) 



(2,0) 



(3,0) 



/ 



3lr 
5 



+ 7j) 



27r 
5 

15 

2 7 2 

?v 

15 " 



\ 



47T 

7_ 
7 



-12949906405.7406 
-12949960422.6225 
-0.00337401832383 



0.0 
0.0 
0.0 



-5651802509.22854 
3767844764.58569 
-0.00147254232049 



-12949906405.7402 
-12949960422.6225 
-0.00337401832383 



0.0 
0.0 
0.0 



-5651802509.22835 
3767844764.58568 
-0.00147254232049 



Note. — The values are evaluated at r = 0.84 R@. The expressions for the 7 coefficients are 
given in Appendix [B] Numerical results are given with 15 significant digits, i.e. up to the numerical 
accuracy. 



where x — (Rt — r) (r — Rb) 2 (R t + — 2r)/{r 2 R^). Rewriting those fields in the conventional spherical harmonics 
writing, we may calculate the axisymmetric components of the vectorial product pXJ x B and obtain 



pU x B 



where the 7 coefficients are defined by: 



27. 
5 

47T 



4tt a/"0 
T Y 3 



37r 

5 



47ryi3 



\ 



2t| 
15 

15 



It 

8?rr 4 
8nr 2 



~5~ 



Te 



l fdsY 3 ° - ( 



47T d e Y? 



(Bl) 



8nR% 



ll 



ll 



and j 2 



9V5/3 5 
8nr 3 R% 

9V5/3 3 R t 



ll 



8nr 5 R% 



■b 1 
1 V - 



These coefficients match exactly the outputs from the code (table 
are then simple scalar products involving the vectorial product 
the analytical calculation. 



if 



The production and flux terms in equation (16) 
They also have been checked by comparison with 



